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In this brief overview we discuss the principal features of real space pairing as expressed via cor- 
responding low-energy (t-J or periodic Anderson-Kondo) effective Hamiltonian, as well as consider 
concrete properties of those unconventional superconductors. We also rise the basic question of 
statistical consistency within the so-called renormalized mean-field theory. In particular, we provide 
the phase diagrams encompassing the stable magnetic and superconducting states. We interpret 
real space pairing as correlated motion of fermion pair coupled by short-range exchange interaction 
of magnitude J comparable to the particle renormalized band energy ~ tx, where x is the carrier 
number per site. We also discuss briefly the difference between the real-space and the paramagnon 
- mediated sources of superconductivity. The paper concentrates both on recent novel results ob- 
tained in our research group, as well as puts the theoretical concepts in a conceptual as well as 
historical perspective. No slave-bosons are required to formulate the present approach. 



PACS numbers: 71.27.+a, 74.70. Tx, , 74.20.-z 

I. INTRODUCTION: APPEARANCE OF 
STRONGLY CORRELATED STATES AND 
PAIRING BY EXCHANGE INTERACTIONS 

The analysis of the unconventional (non-BCS) super- 
conductivity is essentially limited to that appearing in 
the correlated fermion systems. We term the system cor- 
related if the interaction magnitude between the parti- 
cles V is comparable or larger to their single-particle (ki- 
netic, band) energy Eb- In the extreme situation when 
\Eb\ ^ V, we speak about the strongly correlated sys- 
tems. Obviously, this simple theoretical criterion must be 
translated into the specific experimental features distin- 
guishing those systems from other normal metallic, mag- 
netic, and superconducting systems. 

First of all, the short-range strong repulsive Coulomb 
interaction (as compared to the so-called bare single- 
particle energy) will hamper the individual-particle mo- 
tion, and hence increase largely its effective mass to*, 
which in turn, will show up in the strongly enhanced 
value of the linear-specific heat coefficient 7 ^ m* . Also, 
the strong electron-electron interaction leads to the cor- 
responding temperature (T) dependence of the resistivity 
p(r)-p(0) = AT'^, since the coefficient A ~ (to*)^. All of 
these features appear already in the Landau Fermi-liquid 
theory, together with an explanation of enhancement of 
the Pauli susceptibility x ^ ™* ^-i^d in addition, with the 
appearance of collective sound-like excitations. 

The situation has changed decisively with the observa- 
tion of a singular behavior of 7, Xj a-nd A, but with the 
ratio 7/x remaining finite, which appear near the metal- 
insulator (Mott-Hubbard) transition. The transition is 
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associated closely with localization of correlated carriers 
and is signaled additionally by the transition from the 
Slater- to the Heisenberg-type of antiferromagnetic order- 
ing. In such manner, the divergences on the metallic side 
define the borderline of the metallic-state stability. The 
transition occurs for specific (odd) number of electrons, 
e.g. for a half-filled-band configuration of starting elec- 
trons, and thus those systems can be easily distinguished 
from either the band or the Kondo-type insulators, for 
which the fully-occupied bands have even number of elec- 
trons, separated in each case by a gap from empty (con- 
duction) states. Detailed studies of the Mott-Hubbard 
systems near the metal-insulator transitions were carried 
out in the seventies through the nineties of the XX cen- 
tury ^] . They were subsequently supplemented with de- 
tailed studies of orbital, charge, and stripe orderings in 
various Mott insulators, i.e. the Mott-Hubbard systems 
on the insulating side. 

The Mott-Hubbard metal-insulator transitions are 
usually discontinuous, with isolated critical point 
on the higher-temperature side {the classical critical 
point) and a possible quantum critical point at the 
antiferromagnetic-paramagnetic boundary [T]. In a sim- 
ple modeling of this transitions, they are driven by a 
competition between the single particle energy, as repre- 
sented by band energy Eb < (per particle, effectively 
characterized, by bare bandwidth W of fermions) and 
the repulsive Coulomb energy V > (represented by the 
magnitude U of the intraatomic Coulomb interaction). 
The physics, for given number of electrons n per active 
atomic site, is characterized then by changing U /W ra- 
tio (or effectively, by exerting external pressure which 
reduces the U/W ratio). The situation is different when 
we have an orbitally degenerate system, still with one 
electron per active band per site. For example, for n = 1 
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the corresponding Mott insulator may take the form of 
a ferromagnetic insulator with an antiferromagnetic or- 
bital ordering [5] or the state with orbital-selective metal- 
insulator transition. 

A basic question arises what happens if we vary the 
electron concentration (the band filling) instead of chang- 
ing U /W ratio for given n. This situation is quite distinct 
from that when changing U/W. This is because in the 
situation with partial (non-half) filling the metallic state 
is stable even in the strong-correlation limit U/W ^ 1, if 
only the disorder effects associated with e.g. intentional 
doping, do not induce the carrier localization of holes 
or electrons in a weakly doped Mott insulator [3]. The 
doped Mott insulators and the heavy-fermion systems 
are exactly the systems of that type. While the high-Tc 
cuprates such as La2-xSrxCuOi can be represented as 
doped Mott insulator with concentration x < 0.3 of holes 
per CuO\~ active unit, the cerium heavy-fermion stoi- 
chiometric compounds such as CeAl^ or CeCoIn^ can be 
regarded as almost localized systems with 5 < 0.05 holes 
in nominally Ce^^ Af^ electronic configuration (i.e. the 
Ce valence is 3 -I- (5). In both systems the 3d (for the 
cuprates) and the hybridized 4/ — 5d — 6s (for the cerium 
compounds) strongly correlated electrons are itinerant, 
which may transform to the localized states under a mod- 
erate change of stoichiometry, pressure or applied mag- 
netic field. For both systems we assume that U/W ^ 1 
and a; = 1 — n <C 1. 

One may ask whether the clear borderline for the 
change of behavior, corresponding to metal-insulator 
transition, survives also for the doped systems. In other 
words, whether there is a clear distinction between the 
limit of weakly or moderately correlated fermions from 
one side and the regime of strong correlations for n ^ 1 
from the other, when no metal-insulator transition can 
occur. The general assumption usually made is that such 
a dividing line indeed exists (albeit of a crossover type) 
and this statement represents one of the fundamental hy- 
potheses of the theory of strongly correlated systems, even 
though its properties have not been proved convincingly. 
In the case of the cuprates such a line can be drawn 
on the temperature T-doping x plane and terminate at 
the middle point x = Xci ~ 0.15, of the superconduct- 
ing dome and separates a some sort of Fermi liquid ex- 
tends from that line to its upper-end concentration [21 H] 
X = Xc2 ~ 0.3. In t he case of heavy fermion systems, 
the existence of such line is suggested through the ap- 
pearance of antiferromagnetism or metamagnetism in a 
strong applied magnetic field or through the transition 
to the fluctuating (mixed) valence state as a function of 
pressure or else, to the localized moment regime with the 
increasing temperature. 

In our brief overview we assume that such (crossover) 
line exists and thus limit ourselves to the analysis of nor- 
mal and superconducting states in the limit of strong cor- 
relations. This limiting regime is defined as the one, in 
which the probability of double site (orbital) occupancy, 
d^ = {rii^Ui^'j is vanishingly small (so formally, =0). 



This is also the limit, where we cannot start from the 
Hartree-Fock representation of the electronic states as 
the reference state for the further analysis. Instead, 
we will start from the so-called statistically- consistent 
Gutzwiller or Fukushima mean-field approach, devised in 
our group in the last two years [5]. We regard this ap- 
proach as the first consistent renormalized mean-field ap- 
proach of strongly correlated fermionic systems, here em- 
ployed to the description of superconducting state with 
real space pairing. 

One important point should be raised. Namely, the 
pairing interaction is of real-space type and is driven by 
the kinetic exchange interaction under the name of the t-J 
model for a single narrow-band case and the Anderson- 
Kondo model (or the t-J-V model) with the hybrid pair- 
ing for a two-band for a two-band situation. Both of these 
models were introduced by the present author some time 
ago in the context of magnetism and normal correlated 
states and later extended to the description of super- 
conductivity in those systems [HI E] • A somewhat more 
detailed personal account, of the transformation to the 
corresponding effective Hamiltonians into forms express- 
ing explicitly the superconducting pairing, is provided in 
Appendix A. 

The paper is organized as follows. In the next Section 
we discuss the universal aspects of real space pairing, 
namely the fact it is induced by a correlated motion of 
pair of particles throughout the lattice. In Section III we 
summarize the principal features of the so-called Renor- 
malized Mean-Field Theory (RMFT) together with the 
consistency conditions. Particular emphasis is put on the 
phase diagram encompassing the unconventional super- 
conducting states. In Section IV we overview selected 
physical properties of the superconducting state and its 
coexistence with magnetism. In such IV we put into a 
perspective the whole approach. Appendices A-E are to 
provide details of the approach and to put it on a firmer 
formal grounds. 



II. t-J AND ANDERSON-KONDO LATTICE 
MODELS: REAL SPACE PAIRING 

In this Section we overview briefly the universal char- 
acter of real space pairing, i.e. its applicability to both 
high-Tc" and heavy-fermion systems. By the real space 
pairing we understand a correlated electron-pair mo- 
tion (hopping) of partners which are coupled by rather 
strong exchange interaction (i.e. comparable on the 
scale of the Fermi energy). The simplest model of one- 
band strongly correlated metallic systems is the Hub- 
bard model, whereas the Anderson-lattice model reflects 
the corresponding two-band situation when the starting 
atomic (say 4/) electrons are hybridized with the uncor- 
related band (5d — 6s) conduction electrons. Some time 
ago Oil] both of these models were transformed out by us 
to the form expressing explicitly the strong correlations, 
on the low-energy (thermodynamic) scale, by projecting 
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out in a precise manner the doubly occupied atomic con- 
figurations and in this manner reducing the efi^ect of the 
strong intraatomic Coulom b repulsion by replacing them 
with dynamic efi^ects on the low-energy scale. Below we 
summarize briefly each of the two above models. 



A. t-J model and its extension 

Let us start from the extended Hubbard Hamiltonian 
for a single narrow band of s-type: 

^ = 51 *y"Lai<^ + U^ni^nii + 7;^ K^jn^Uj, (1) 



where the primed summations mean that we take only 
i ^ j terms in both the hopping (the first) and the in- 
tersite Coulomb (the third) terms. The second term is 
the celebrated Hubbard term representing the energy in- 
crease when the site (Wannier) state is doubly occupied. 
The principal feature here is that the Coulomb repulsive 
interactions are developed systematically and express re- 
spectively the interaction between two electrons on the 
same site (with magnitude U) which represents (by far) 
the largest energy scale, as well as between the electrons 
on neighboring sites i j, with magnitude Kij and ne- 
glect all more distant, interactions, since we assume that 
Wannier states are strongly localized ( the tight-binding 
approximation). In other words, this atomic representa- 
tion differs drastically from the concept of a lattice elec- 
tron gas as a starting point, since the long-range nature of 
the repulsive Coulomb interaction is cut off. The decisive 
step was taken [S] to derive an effective Hamiltonian out 
of expressing low-energy dynamics, i.e. which con- 
tains high-energy processes as virtual transitions in the 
second order order. The resultant effective Hamiltonian 
in the lowest Hubbard band is: 
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The double-primed summation means i =/= j =/= k =/= i. 
The projected fermion operators are: b]^ = al^ (1 — Uig), 
bja = fljo- (1 - Jijff), I'lG = b\^bia, Vi = I^^i^icr, and 
5, = (5f,5,f) = {b\^h-,\v,^-ViC)l'^- We see that 
apart from a restricted (projected) hopping, with no dou- 
ble occupancies present, we have also present an antifer- 
romagnetic kinetic exchange interaction with exchange 



integral 



4if,/(C^ - Kij) [8] (the second term). 



here generalized to the case of itinerant fermions, since 
the spin operators are explicitly expressed in terms of 
fermionic operators. The third term describes three-site 
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FIG. 1. (Color online) Top: schematic representation of the 
Hubbard subbands. The virtual hopping processes corre- 
spond to the virtual transitions lower — >■ upper — >■ lower Hub- 
bard subbands (process a). Bottom: three possible intersite 
hopping processes (in real space) determining the dynamics 
of Hamiltonian (2) in the lower Hubbard subband (i.e. for 
the band filling n < 1): a) virtual hopping between single oc- 
cupied neighboring sites that leads to their antiferromagnetic 
(kinetic) exchange interaction; b) single-particle hopping be- 
tween the single occupied and the empty sites; and c) three- 
site hopping between the singly occupied and the empty site 
via an intermediate single occupied site with opposite spin. 



hopping processes, and the last - the residual repulsive 
intersite Coulomb interaction. The dynamical processes 
in the second order, taken into account in ([2]) are shown 
in Fig 1. We also included there the renormalized single- 
particle hopping. 

One has to mention that the projected fermion op- 
erators {^io-} and do not obey the usual fermion 
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anticommutation relation, as we have that 

{ha.b]^,} = S,^ [(1 - n,^) 5,,, + Sf (1 - S,^,)] 



(3) 



This last property creates a basic formal complication, 
as we are working now with projected fermion operators. 
We can say that the effective model represents, in prin- 
ciple a fermionic quantum liquid, which is not a Landau 
Fermi liquid (usually termed imprecisely as: non-Fermi 
liquid). As a first approximation, we discuss in the next 
Section an effective {almost localized) Fermi liquid and 
demonstrate the usefulness of the approach as a renor- 
malized mean field theory. Note that the Hamiltonian 
([2|, is valid for the band filling n < 1. It can be easily 
extended to the situation when n > 1 by employing the 
hole language. 

The model ([2]), mainly in the limit of Kij = 0, was 
used to discuss the magnetic, charge- density, and mixed 
(stripe) phases. However, a new impetus to study the 
model was the idea [3], that the t-J model (i.e. Hamil- 
tonian ^ with Kij = 0) can be used to describe the 
so-callea reaZ space pairing. The real-space pairing oper- 
ators may be defined in the following manner 



4- V2 



^ (b^ 



Bij = ^ {bi^bj^ - bi^bji) . 



(4) 



In this manner, both two- and three-site terms can be 
recast to the closed form and the effective Hamiltonian 
transforms to the form 



2i? 



ij <T ij k ^■^ 



(5) 



where now both the terms with k = i and k ^ i are 
incorporated into the second term. The expressions Q 
represent the projected spin-singlet creation and annihi- 
lation operators located on the pair of sites {i,j) with 



i ^ j, since B, 



Bi. 



0. Hence, no ionic mixture to 



such two-site spin-singlet state appears. Also: 
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(6) 



i.e. the "number operator for local singlets" is equiva- 
lent to the generalized Dirac exchange operator (r.h.s). 
In other terms, the description in terms of itinerant-spin 
interaction language is equivalent to the description in 
terms of local {i,j) itinerant-spin singlet. In the other 
words, the antiferromagnetic ordering coming from the 
kinetic exchange interaction and the kinetic-exchange or 
real-space singlet pairing should be regarded as equiva- 
lent ways of describing the interaction i.e., express a dif- 
ferent type of ordering, which should be treated on equal 
footing. What is even more important, such paired state 



is directly included by the correlated motion of those sin- 
glet pairs, as discussed below. 

One general remark is in place here. The operators 
Bjj and Bij express explicitly the singlet nature of the 
pairs if {Bjj) ^ 0. Normally, one operates only with 
either (cl^Cj^) or alternatively, with (cj^c^^) [TBI. But 

then, one implicitly assumes that (cLcj^) = 
Strictly speaking, this relation should be checked out ex- 
plicitly, e.g. evaluating those two averages separately. 
Nota bene, this explicit checkout helps in distinguishing 
between the singlet pairing and the triplet pairing (both 
components with = 0). Namely, in the latter case 
(cl^Cj^) — ^P™ nature of real-space pair- 

ing is particularly obscured, if we consider the coexistence 
between magnetism and superconductivity, as discussed 
below. 

Additional question is concerned with the presence of 
the repulsive term ~ K^j, which is usually neglected 
in the analysis of high-temperature superconductivity 
within the t-J model. Omission of this term is justified 
by the circumstance that is regarded as contributing only 
the reference energy, i.e. is the same for the phases under 
consideration. It may be also regarded as compensated 
by the electron-lattice interaction which leads effectively 
to Kij < [17]. If this would be the case, then that the 
pairing part induced by the case Jij does not lead to iso- 
tope effect, whereas the part ~ Kij < does. Since the 
isotope effect in high- Tc system is small, it means that 
the kinetic exchange part of pairing in dominant. 



B. Hybrid (Kondo-type) pairing in 
Anderson-lattice model 



The Anderson-lattice model represents the simplest 
two-orbital model with a coherent mixture (hybridiza- 
tion) of strongly correlated and uncorrelated electrons. 
It is the canonical model of hcavy-fermion system if the 
orbital degeneracy of e.g. 4/^ electronic configuration of 
Ce'^+ ions is not crucial (e.g. the lowest double is the 
only important crystal-field level which hybridizes with 
conduction electrons originating from the 5d — 6s itiner- 
ant states). The starting Hamiltonian has the following 
form in the site (Wannier) representation: 

inner ia i 

+ E (^™aLcmtr + Vi^C^^^aia^ + Ufc ^ N^'^rn, (7) 

where (i,j) label starting atomic (a) states, (m,n) label 
starting delocalized (conduction, c) states, and = 
^\(y^i<y represents the hybridization matrix ele- 

ment. The strongly-correlated aspect shows up through 
the circumstance that both \Vim\ and \tmn\ ^ U, but 
the position of the atomic level £/ ~ 1^™, so not the 
whole hybridization term can be transformed out via the 
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SchriefFer- Wolff-type of transition to the Kondo lattice 
type of model [71 [T9]. 

To adopt the model to the situation with strong corre- 
lations, one performs the transformation, in which as in 
the narrow-band case, the double occupancies of atomic 
states are excluded and replaced by the virtual processes 
leading among others to the Kondo interactions between 
the electrons of the two subsystems J_^. Explicitly, the 
full effective Hamiltonian in the second order in V/U and 
representing the low-energy dynamics, takes the form: 



U + ef 



+ H.C 



.... ^ + 'f 



(8) 



where 

|'„<T Z 



aLa. 



ml; 



(i-iv^a), 

?^m;)/2). 



and Sm = 



The important point to note is that in this Hamil- 
tonian we have both antiferromagnetic Kondo coupling 
with the exchange constant Jim = 2\Vim\^/{U -t- e/) 
and the residual hybridization in the projected subspace, 



-\- H.C.). In this manner, the itineracy of 



/ electrons is explicitly expressed unless there is a phase 
transition to the localized Mott state for / electrons (i.e. 
to the state with (i/^) = 1, see below). The dynamic 
(virtual) processes in the second order and taken into 
account in the effective Hamiltonian are displayed 
in Fig. 2. There, we also show explicitly the remaining 
processes: hoping in the conduction band and the ef- 
fective /-electron hopping (the bottom most part). The 
second important point is, that in analogy to the real- 
space spin-singlet pair operators ([5]), we can introduce 
the hybrid (Kondo-type) pair operators, namely define: 



V2 ( 



, J — J 

^ (fe^Cmt — bi^Cml) ■ 



(9) 



and rewrite the Hamiltonian in a more closed form, at 
least when Vim = ^imi neglect the second order 

contribution to the hopping (the second part in the first 
term of ([7|). We obtain then 



(10) 



+ XI (Ccma + cl-,J)^^ - X ^"^'J' bXmb 

ima imn •' 

Note that, in analogy to ^ we have now 
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FIG. 2. (Color online) Top: schematic representation of 
hybridization-induced process as f-level occupation dependent 
hopping process and its division into the low- and the high- 
energy processes. The former (I) corresponds to the presence 
of the residual hybridization term in (8); the other (II) leads 
to the Kondo-type coupling which in turn is expressed as a 
real-space hybrid pairing in (11) in the second order in V/e/. 
Bottom: the possible interband hopping process in real space 
involving virtual hopping process to the double occupied /- 
level a) and the three-state (site) hopping process ( b) and c) , 
2 and 2' alternative steps). The arrow (labels) by tij repre- 
sents the resulting intersite f — f hopping process. 
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i.e. both the Kondo interaction and the hybrid singlet 
pairing in real space appear on the same footing. Both 
processes are characterized by the Kondo exchange in- 
tegral J4 = 2\Vi.m?/{U + e/). This is the Anderson- 
Kondo Hamiltonian capable of describing both the itiner- 
ant nature of the heavy-/ electrons and their localization, 
as well as the hybrid (Kondo) paired state in the limit 
of strong correlations among the hybridized fermions. 
By analogy with t-J model (5]), the one represented by 
Hamiltonian (10 1 can be cal ed the t-J-V model. Also, 



when going to the higher, fourth order in Vim/iU + £/), 
we can include additionally the /-/ pairing via the term 
~ V^/{U + effB\.Bkj. In effect, both the (orbital se- 
lective) localization of / electrons, as well as a two gap 
superconducting state are encompassed as limiting cases 
of this coherent quantum liquid composed of two sets of 
hybridized fermions. 

To summarize this Section, the formal expressions 
of the kinetic exchange and of the Kondo interaction 
through the same type of real space pairing operators 
illustrates the universality of the pairing in both the 
Mott-Hubbard and heavy-fermion systems. Obviously, 
the practical usefulness of the pairing concept is proved 
only by showing the stability of the corresponding mag- 
netic and superconducting states induced by this uncon- 
ventional forms of pairing. In the next Section we sum- 
marize formulated by us [HllS] the statistically-consistent 
rcnormalized mean-field theory, as well as present some 
exemplary results. 



III. RENORMALIZED MEAN-FIELD THEORY: 
STATISTICALLY CONSISTENT APPROACH 
(RMFT-SCA) 

A. Gutzwiller approximation with 
statistical-consistency conditions(SGA) for the 
Hubbard model 

The above Hamiltonians have a complex form, since 
they contain the projected fermionic operators (i.e. com- 
posite fermion operators with non-fermion anticommu- 
tation relations). The question is how to diagonalize 
Hamiltonian containing such operators, at least in an ap- 
proximate and consistent manner, so one has the confi- 
dence of having a well defined mean-field approach. Par- 
enthetically, note that by making the canonical trans- 
formation we include a certain class of higher-order dy- 
namical processes automatically. Hence, a relatively sim- 
ple approximation (of the Hartree-Fock type) on the 
effective-Hamiltonian level includes, at least partially, 
those higher-order dynamic processes. Therefore, instead 
of rigorous solution impossible to attain for many-particle 
correlated systems with spontaneous symmetry break- 
down, we develop consistency check for approximate so- 
lution at hand. 

Actually, we can do a bit better than just carrying out 
the Hartree-Fock type decoupling in the effective Hamil- 



tonian. Namely, the name of the game is the renormal- 
ization of the effective Hamiltonian combined with a sub- 
sequent Hartree-Fock decoupling of the many-body parts 
still remaining. This procedure is not systematic in the 
same sense, as the canonical perturbation expansion pre- 
sented in the preceding Section, but it is in our view , in 
accordance with our physical intuitive insights into the 
nature of strong correlations. 

Probably, the simplest nontrivial approach to the de- 
scription of a correlated state is the Gutzwiller approach, 
discussed here in its still simpler form of a Gutzwiller 
ansatz. It relies on a variational approach by postulating 
the corresponding macroscopic (A'^-body) wave function 
\^) and a subsequent approximate combinational evalu- 
ation of the relevant averages appearing in and 
in (5" 1 5"). The wave function respects the exclusion of 
the double occupancies in real space, in an approximate 
manner though. An altered approximation scheme has 
been introduced subsequently by Fukushima [10. who in- 
troduced additional variational parameters, the so-called 
fugacity factors, which guarantee that expectation val- 
ues of the number of particles in the uncorrelated and 
correlated states are equated as the same. 

In the Gutzwiller approximation the postulated wave 
function has the form 



(12) 



where [vfo) represents the Fermi sea of uncorrelated (usu- 
ally noninteracting) electrons, and <? is a variational pa- 
rameter equal to zero when the double occupancies are 
excluded. By applying this type of wave function to 
the narrow-band Hamiltonian, we can obtain the ground 
state energy Eq/N per site of the spin polarized state in 
the form 



Eg 
N 



^zZ^'^ (d, ria) ecr + Ud^ 



(13) 



where d^ = {ui-fnii), e^ is the average band energy of 
particles with spin a (per site), n„ is the corresponding 
average number of particles with spin u, and 



(la = 



{[{n, -d^){l-n + d^)] + [d^ {n, - d')] 



(1 - ria 



(14) 

is the so-called band narrowing factor which represents a 
renormalization factor, ^ ^ 1, of bare band energy 
under the influence of correlations (in the Hartree-Fock 
limit — {ni„){nia) or equivalently, = 1). Alterna- 
tively, one can write down ( |13| as an expectation value 
for the bare wave function of the following effective 
single-particle Hamiltonian [TT]: 

"Hga = ^ qa{n, d, m)tijc\^Cja + NUd^, (15) 



where w = an^ is the spin magnetic moment per 
site. This effective Hamiltonian contains two terms: the 



7 



renormalized hopping and the expectation value of the 
Hubbard term. The replacement of the Hubbard term by 
its expectation value means that we consider the single- 
particle propagation in a frozen configuration of the dou- 
ble occupied sites, which at the end is optimized by min- 
imizing Eg with respect to d. 



the corresponding minimum condition for J^, what guar- 
antees automatically the statistical consistency as one 
can see by imposing dJ^/dXm — dF /dX^ = 0. 

One should note that the variational parameter d is 
not a thermodynamic variable, so the corresponding con- 
straint does not appear. Nonetheless, in such modified 



Hamiltonian (15 1 can thus be easily diagonalized and formulation all the variables appearing in (16) are calcu- 



the corresponding free-energy functional J- constructed 
in the standard manner has the form: 



j,(GA) 



ga^y^ In 

ka 



(16) 

The quantities m and fi can be calculated either by min- 
imizing this functional for selected n, T, and U/W or al- 
ternatively, by writing down selfconsistent equations for 
them, i.e. 

iY.^'^kaaka); ^Y.'^i'^L'^ka)- (17) 



Obviously, (al^cika) = "-fccr is the average occupancy - 
the Fermi function. 

In the executing of either of the two above procedures 
one can show explicitly that they do not yield the same 
results. This means that the fundamental principle of 
Bogoliubov, which holds for the Hartree-Fock approxi- 
mation, is not obeyed here. This problem can be traced 
to the presence of renormalization factor Qo- in either Eq 

or J", since then Y.ka gi'm"''"'^ <^k{nka) ^ 0. This very 
important difficulty has been omitted in almost all pa- 
pers utilizing the Gutzwiller approximation. 

To overcome this difficulty, we have proposed [5] that 
the renormalized Hamiltonian ( 15 ) has to be supple- 



mented with additional constraints, expressed with the 
help of Lagrange-multiplier method when calculating the 
averages of the type (17), providing an appropriate con- 
dition minimum of thus corrected Landau functional J- . 
Note that such modification will concern any renormal- 
ized in this manner mean-field Hamiltonian (e.g. the 
above t-J or Andcrson-Kondo models) as it contains 
nonanalytic renormalization factor q„. In this manner, 
the statistical consistency of the whole approach (in this 
sense of a correct statistical physics) is guaranteed, i.e., 
the self- consistent equations and the variational mini- 
mization provide the same results. Explicitly, in the sit- 
uation when we consider normal state, i.e., when only m 
and /X (or m and n) appear as a thermodynamic variables, 
we define now the effective Hamiltonian as fo Hows 

a = 'HgA - I ^ crnfc^ - ^ (T{nka) j 
\ k(7 ka / 

\ ka ka / 



(18) 



where the Lagrange multipliers and A„ play the role of 
extra (global) molecular fields to be calculated also from 



lated variationally and thus the approach becomes self- 
consistent and self-contained. For the sake of simplicity 
we do not elaborate in detail of the approach using the 
Fukushima approach |13H15| . However, we discuss some 
of the results obtained by using this particular method. 



B. Statistical consistency for t-J model and 
Fukushima variational wave function (SCA) 

The statistical consistency conditions have been im- 
plemented practically simultaneously to the two models 
interesting us in the context of unconventional supercon- 
ductivity: the extended t-J model with real space pair- 
ing induced by the antiferromagnetic kinetic exchange 
[SI [131 [2] as well as to the Anderson-Kondo model with 
hybrid real space pairing induced by the Kondo coupling 
[TSl [T^ . Both of these types of pairing have been intro- 
duced in the preceding Section. 

The situation in the t-J model has been additionally 
modified by taking into account a modified Gutzwiller- 
type wave function introduced by Fukushima IH] . In the 
latter approach, the Gutzwiller wave function ( 12 ) in the 
strong-correlation limit {g — 0) is replaced by 



I*f) = n iMT'''^ {HT"'^ (1 - «^T"4) l^o) = Pf\^^) 

(19) 

Such modification of the Gutzwiller projector Pg allows 
us imposing the condition, that uncorrelated average 
number of electrons per site (on each site "i") and with 
spin is equal the corresponding actual average computed 
within the scheme. In other words, 



(20) 



where subscript " F " mean evaluation with the wave 
function (|9|. This conditions represent additional self- 
consistency requirement. In effect, as well will see, some 
of the renormalization factors are the same, some are dif- 
ferent. 

Two important practical points concerning this 
method should be emphasized. First, since the ground 
state energy is expressed as 



Er 



(21) 



and l^o) represents an uncorrelated state, the multiple- 
fermion-operator expectation values contained in ( 20 ) 



can be evaluated by factorizing them into those contain- 
ing only pairs of operators, each evaluated for uncorre- 
lated state (a Wick-type contraction!). This is somewhat 
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cumbersome procedure, so we present only the final re- 
sult: 



A,: 



1 - ("-ia) 



1 - (rii) ' 

as well as the form of the effective t-J Hamiltonian: 
(^-oIPfHPfI^o) 



(22) 



n 



- E4''^[(c-Cj--^»J-) + fi'-c.], (23) 

where the effective t-J Hamiltonian SCA H is expressed 
by a corresponding expression SCA, which will not be 
reproduced in detail here (see [13]). However, on should 
note that in ( 23 1 we take in this method an expectation 



value of (H). Thus the operator part of the Hamiltonian 

is composed of the constraints. Such trick is analogous 
to that leading to the effective single-particle Hamilto- 
nian (17), where the expectation value of the interaction 



part is taken. Note also, that we evaluate {Ppl-LPp); as 
well as {Pp) as expectation for the noninteracting state 
|\l/o). So, we can utilize the Wick-type factorization of 
the operators in the real space representation. 



C. Statistically consistent approximation for the 
Anderson-Kondo lattice model 



Very recently, we have extended [191 the Fukushima 
approach to the periodic Anderson-Kondo Hamiltonian 
(10). First, we have extended our former results [20] for 



heavy-fermion systems to the case with nonzero applied 
magnetic field. Second, the work is progressing on the 
hybrid real-space pairing and its coexistence with mag- 
netism. The latter work |21j extends our earlier work on 
unconventional superconductivity with hybrid pairing in 
a lower hybridized band and in the Gutzwiller approxi- 
mation. This work will be elaborated elsewhere [TSj |19J . 
Also, it would be desirable to analyze high-temperature 
superconductivity within d—p hybridization included ex- 
plicitly and within the SCA approach. 



D. Statistically consistent approximation: A brief 
outlook 



It is very important to note that the constraints intro- 
duced and composing an essence of the SCA scheme are 
the same, at least for the normal state, with the more 
involved auxiliary (slave)-boson type of approach in the 



saddle-point approximation [5]. However, our approach 
contains no auxiliary ("ghost") Bose condensed fields, 
which introduce spurious phase transitions. Second, the 
present approach is quite natural on physical grounds, as 
well as extends the long tradition of introducing a molec- 
ular field as a conjugate variable to each introduced order 
parameter, here in the situation with a complex ordering. 
Also, as one can see already from the complex form of the 
mean-field Hamiltonian, that it represents a quantum liq- 
uid state intermediate between the Landau Fermi-liquid 
state and what is called a non-Fermi liquid state, in ad- 
dition to providing the Mott or selective-orbital-Mott lo- 
calized states, as quite unconventional, as we discuss it 
elsewhere. 

In summary, the SCA incorporates the slave-boson 
approach into a statistically-consistent single-particle 
mean-field theory with a complex formal structure in- 
volving a number of order parameters and associated 
with them mean fields. Once the form of a variational 
wave function is selected, all the remaining analysis can 
be systematic in the sense of many-particle perturbation 
theory [22j. 



IV. PHYSICAL PROPERTIES: 
SUPERCONDUCTING AND MAGNETIC 
STATES 

The most important feature of any theory is to provide 
a quantitative, or at least a coherent semiquantitative, 
description of the relevant electron states and properties 
and in particular, a phase diagram involving physically 
plausible phases that given fermionic model allows for, at 
least in the in mean-field approximation. Below we char- 
acterize separately recent results for the t-J model and 
for the hybridized (two-orbital) Anderson-Kondo model, 
in both cases involving the stable phases. 



A. High-temperature, single plane 
superconductivity within t-J model in SCA 
approximation 

To determine the role of different terms we rewrite the 
Hamiltonian ^ in the following form 



ija (ij) 



(24) 



where ci 



C2 



1 and K 



AtjAU — Kij), represents the three-site terms, and 
{ij) means the summation over nearest neighboring pairs, 
each taken once. The factor ci and C2 have been intro- 
duced to differentiate between the various contributions 
to the paired state energy. Also, we consider here only the 
case with Kij = 0, i.e. neglect the intersite part of the 
repulsive Coulomb interaction. The rationale behind the 
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last assumption is that, since we consider a correlated- 
pair motion as the source of pairing, the repulsive contri- 
butions is only one of the contribution, not the decisive 
one. Additionally, we discuss here only the final results, 
as the formal analysis has been elaborated in detail else- 
where [T3UI5] . Let us mention only that the result can 
be represented, by the Hamiltonian (24) reduced to the 
renormalized form 



(ren) 
t-J 



ijij 



A*-A- ■ 



where 



(25) 

glj and ^0 are the corresponding factors. The form 



of follows from a detailed procedure which leads to the 
renormalized hopping in the correlated (I^'f)) state: 




(26) 



where Xiju = i'Ha'^jcr) ^^'^ hopping amplitude in the 
uncorrelated state and therefore for the paramagnetic 
state with real gap parameter Ay- 
have that 



= (4cj-) = A*^. we 



9h 



and 



1 - n. 



1 — rii 



Xi 



A2. 



A 

9^3 



1 - Tli 



1 — rii 



A2. 



(If) (1 



(?7) 



.(28) 



Here rii = (rii). These renormalization factors reduce to 
the standard Gutzwiller-ansatz form for spatially homo- 
geneous case, rii = hj = n. Here, to first approximation 
9ij — 9iji the formation of the pairing requires elec- 
tron itineracy. Note, that strictly speaking 



(ren) 



-J 



(29) 



Therefore, the exchange and other interaction parts not 
written down explicitly, should be decoupled in the mean 
field manner using all possible (assumed as nonzero) con- 
tracted pair-operator averages. In effect, we can diago- 
nalize thus approximated Hamiltonian ( |24[ ) with the de- 
coupled terms and supplemented with the constraints. 
The resultant renormalized hopping amplitude and 
the gap parameters Dt in the r-th direction of the lat- 
tice, have the following forms: 



Tr=trgl + -Jr9rXr+\'^^\ 



(30) 
(31) 



where we have included only the terms written down ex- 
plicitly in (25 1. Also, we have assumed that Xx — Xy — X 



,Sr,CuO, 



non-Fermi liquid 




0.02 0.06 0.2 0.32 

doping level (holes per CuOj) 
underdoped optimally doped overdoped 



0.035 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 
DOPING, X = 1-n 

FIG. 3. (Color online) Top: overall phase diagram for 
La2-xSrxCuo4 on the temperature T - doping x plane. Bot- 
tom: Doping dependence of the renormahzed superconduct- 
ing order parameter Ac. The curves 1-7 are explained in de- 
tail in Ref. [11] and correspond to various methods of solving 
Hamiltonian (2). The most important feature of the solution 
is the appearance of the upper critical concentration for the 
d-wave superconductivity disappearance. The vertical line 
defines roughly the optimal doping. 



and Aj; = — Aj, = A, the second relation being for the 
superconducting gap of the d-wave form. 

In Fig. 3a we plot the overall phase diagram for the 
high temperature superconductor La2~xSrxCui. The 
corresponding theoretical phase diagram containing the 
superconducting part of the phase diagram is exhib- 
ited in Fig. 3b. Few features of theoretical phase dia- 
gram should be noted. First, the superconducting dis- 
appears only in the Mott insulating limit n — >■ 1 (dop- 
ing x = 1 — n — 0). This inadequacy of our approach 
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are caused by the circumstance, that we have not in- 
cluded in our analysis either the onset of antiferromag- 
netism or the atomic disorder induced by the Sr dop- 
ing. Second problem with our approach is connected 
with the absence of the pseudogap appearance. This is 
because we neglect here the phase fluctuations of the or- 
der parameter Aij. Nonetheless, the approach contains 
the two very attractive features. First, we show clearly, 
that there exist an upper critical concentration for the 
disappearance of superconductivity. The various curves, 
which were calculated with different approximations in 
the Hamiltonian (as discussed in detail in [T3]), show 
clearly that the upper critical concentration is the dop- 
ing regime x — 0.22 — 0.35, in agreement with the ex- 
perimental data for various single-plane cuprates. One 
should state right away, that the existence of this critical 
concentration speaks out decisively in the favor of real 
space pairing, as with increase of doping the hopping in- 
creases, which in conjunction with dilution electrons with 
X destroys the correlated pair motion. Second, here we 
have a well defined Fermi-liquid state, particularly in the 
so-called overdoped regime (to the right of the vertical 
line in Fig. 3b). The real space involves a correlated-pair 
motion and has been animated by us recently [2B] . To 
put it bluntly, the pairing discussed here does not involve 
any intermediate boson such as paramagnon. Finally, the 
magnitude of the renormalized gap Ac = g^{Aij) is in 
the unit of the first hopping integral \ti \ which is taken in 
the range 0.35 — OAOeV {\ti\/J = 3 is assumed). There- 
fore, the maximal value of A ~ 100 — 120K is obtained 
for the uppermost curve, a quite reasonable value in view 
of the mean-field nature of our estimate. An elementary 
interpretation of the superconducting "dome" shape de- 
picted in Fig. 3 is based on the circumstance that here 
renormalized hopping amplitude ^ \t\x is comparable to 

J ESI. 
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FIG. 4. (Color online) Doping dependencies of the SC gap Dk 
at it = (tt, 0) for cases 1-6 and for t' /t = -0.27, and J/\t\ = 0.3 
(filled diamonds). Large filled circles - experimental data. For 
a detail discussion see Ref. |13) . 

In Figs. 4 and 5 we display respectively, the doping 
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FIG. 5. (Color online) Doping dependence of Fermi velocity 
in the nodal ((0, 0) — >■ (tt, tt)) direction. For a detail discussion 
see Ref. [13]. 



dependence of the quasiparticle energy in the antinodal 
direction (fc = [tt, 0]) and the Fermi velocity in the nodal 
([0,0] — >■ [tTjTt]) direction. The points represent various 
experimental results for the single-planar systems |13| . 
Our analysis represents a more precise version of that 
presented in Ref. |14j . 

From Figs. 3-5 we can see that our results express in a 
semiquantitative manner the trend of the corresponding 
experimental data for overdoped systems (to the right of 
the vertical line). However, it must be reiterated again, 
that the phase diagram does not contain the crossover 
(pseudogap) line, cf. Fig. 3b, as well as provide a wrong 
tendency of the Fermi velocity wf, i-c, its almost con- 
stant value in the underdoped regime. This most proba- 
bly means that pocket of Fermi-surface (the arcs P5| ) 
near the nodal direction survive almost untouched in 
underdoped regime. The situation in the underdoped 
regime requires a generalization of the present SCA ap- 
proach (some hints to that generalization are provided in 
Appendix E). 

One more interesting feature of our results should be 
noted. Namely, the temperature dependence of the su- 
perconducting gap magnitude Ac = Ac(T) follows to 
a very good approximation the BCS dependence in the 
version proposed by Rickayzen |26| , which in our context 



Ac(r) 

Ac(0) 



tanh f ^ 
VA(0) 



T 



(32) 



with relative temperature t = T/Tc- The comparison 
our formula (32) is illustrated numerically in Fig. 6. This 



is an amazing result given the complicated nature of the 
self-consistent equations leading to the numerical results. 
For the sake of completeness we derive the above equation 
in Appendix B. We also show there that 2/S.c{^) I^bTc — 
4. 
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FIG. 6. (Color online) Reduced temperature dependence 
of renormalized superconducting order parameter {A)c = 
Ac(T') for various doping levels: squares: x — 0.125, solid 
circles: x = 0.25, diamonds: x = 0.3 are all for ci = 1, 
C2 = 0, and t' /t = -0.25. Solid line: BCS result ((32 I. 



B. t-J model: coexistence of magnetism and 
superconductivity 

A complete analysis of the t-J model requires determi- 
nation of the Fermi-surface evolution with doping within 
this almost localized Fermi-liquid picture, along the lines 
presented in |13j . Furthermore, one should incorporate 
the magnetic phases into the phase diagram presented in 
Fig. 3b. The last type of analysis has been performed 
very recently [HI [18] and the results are presented in 
Fig. 7, where the same type of approach as above has 
been applied with one addition feature. Namely, a mix- 
ture of the spin-singlet and the so-called staggered spin 
triplet states appears when the antiferromagnetism co- 
exists with superconductivity j27| . The appearance of 
the triplet component of the gap is imminent in that 
situation even though we have only the singlet-pairing 
interaction or equivalently, when only the antiferromag- 
netic exchange interaction in present. This is because of 
the following reason, illustrated in Fig. 7. Namely, if we 
have two-sublattice antiferromagnet, the pairing ampli- 
tude Aa = (ci^Cji) between the spin majority electrons 
sublattices {i S A,j S B) should be different (larger) 
than that on the corresponding quantity = (ci^Cj^) 
for the minority-spin electrons. Additionally, we would 
like to describe the d-wave superconductivity, expressed 
in the real-space language, which amounts to postulating 
the following form of the gaps 



TijAA for ie A 

t„ Ab for ^ e S ' ^ ' 



where Tij = +1 for the n.n. pair (ij) along x and Tij = —1 
along y axis, respectively. One can note immediately, 
that if (cjiCi-f) and (cj-fCn) averages differ in magnitudes, 
then we can have a mixture of singlet and triplet pairing, 
as in e.g. singlet-pairing case they would be of opposite 




A 



V 



FIG. 7. (Color online) Spin-majority (blue, bigger arrows) 
and spin-minority (red, smaller arrows) electron spins in a 
system with the AF order and real-space superconducting 
gaps. Aa binds two spin- majority electrons, and A_g binds 
two spin-minority electrons and therefore, there is a priori no 
reason for these two gaps to coincide (as would be the case for 
no staggered yr-triplet component). In other words, the two 
distinct gaps make effectively the t ~ 4- ^nd i — t pairing 
components of the opposite-spin pairs distinguishable. 



sign and for pure triplet the sign would be the same. In 
effect, the unrenormalized gap parameter A^j can in the 

present situation be decomposed into the singlet (A^p) 

and the triplet (A^ ) parts according to a prescription 



A ■ = A^'^^-l-A'^'e*^-^' 



(34) 



with 



A 



(s) 



Af 



AT) 



Af 



(35) 

where Q = (tt, tt) is the superlattice vector in the case 
of simple two-sublattice AF ordering in two dimensions. 
The amazing feature of the present representation of 
the gap in the form ([35]) is that the triplet component 
is present even though there is no explicit spin-triplet- 
pairing inducing interaction in the Hamiltonian. The 
detailed analysis has been performed in |14j . where the 
corresponding statistical-consistency constraints and the 
minimization of the (appropriate free-energy) functional 
has been carried out. One should note that in the present 
situation we encounter, for the first time in this approach, 
the spatially modulated occupancy, i.e, it has the form: 



(^ol^i^l^o) 



I: {n + amp + (JmAMe"^'^") , 



(36) 

where, as before, n is the band filling (0 ^ n ^ 1), mpM 
is a ferromagnetic (homogeneous) spin-moment compo- 
nent, and mAF is the antiferromagnetic (staggered, sub- 
lattice) moment component. The results, encompassing 
a simultaneous optimization of a system of 11 algebraic 
equations for the case of square-lattice case, is shown in 
Fig. 8 in the form of the phase diagram on the applied 
magnetic field h (in units of \t\) - band filling plane. The 
microscopic parameters taken in that computation are 
listed in that figure. 

Particularly interesting for us is the horizontal line, 
where we observe the sequence of phases AF + SC 
AF — PP — FM with the decreasing band filing. In a 
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Phase Diagram: T = 1/|3 = 0.002, t'=t/4, J=t/3=1 




0.98 0.96 0.94 0.92 0.9 0.88 0.86 0.84 0.82 0.8 
Band filling, n 



FIG. 8. (Color online). Phase diagram on the band filling 
- magnetic field plane. The phases are labeled as follows; 
AF+SC - phase with coexisting superconductivity and anti- 
ferromagnetism, AF - antiferromagnetic phase, FM - ferro- 
magnetic phase, SFM - saturated ferromagnetic phase (with 
rriFM ~ n). For further analysis we restrict ourselves to 
n = 0.97 as marked by the dashed vertical line. No stable 
pure superconducting solution has been found. 



way, the results are somewhat disappointing, as we would 
rather expect, with the increasing n, a clear transition 
from AF Mott insulating state to a pure superconducting 
state. Here the coexisting phases AF + SC phase are 
stable up to Xc' = 0.0.6, and the for x > Xc the pure 
two-sublattice AF state becomes stable. 



C. Fulde-Ferrell phase in narrow-band limit 

In this Section we overview briefly plausibility of the 
observations of the Fulde-Ferrell phase in the strongly 
correlated-electrons systems, as viewed from the point of 
view of our SG A approach f28^ . 

The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) was 
proposed theoretically many years ago [59]. The mo- 
tivation for our work in this area was the suggestion 
that this state has been observed in the heavy-fermion 
system CeCoIn^, possibly coexisting with antiferromag- 
netism [5CTj . It was quite a coincident with the experimen- 
tal observation of spin-dependent heavy electron masses 
|31j . On one hand, the appearance of the spin-dependent 
heavy electron mass was proposed earlier as one of the 
crucial phenomena for the strong correlated systems [35] . 
Therefore, a natural idea appeared to discuss the effect 
of the spin-dependent masses (SDM) on the FFLO state 
stability. However, we discuss first briefly the concept of 
the spin-dependent masses. 



Spin-dependent masses 

For a quasiparticle gas of correlated particles, their en- 
ergy in the applied field h = ^gfisHa, when the Landau 
quantization is neglected, is expressed as (2S] 

Ckcr = ah- ^1- ahcorr, (37) 

2m.a 

where h^orr is the effective field induced by the correla- 
tions and spin-dependent mass enhancement ma-fmB is 
of the form in the limit of C/ — >■ oo for a single narrow 
band [31] 

TTicr 1 — rio- 1 — n/2 m 
TUB I — n 1 — n 2(1 — n) 

= — {ruav - ^Am) , (38) 
mriB ^ 2 / 

where tub is the bare (band) mass and to = rt^ — is 
here the system spin polarization. Also, Am = to^ — 
is the mass difference, while niav = {tti^ + m\) /2 is the 
average quasiparticle mass, i.e., the mass in the absence 
of magnetic polarization. It is interesting to note that in 
the magnetic-saturation limit we recover the band limit 
for the spin-majority subband, i.e. to^/tob — 1, whereas 
the heavy quasiparticles in the spin-minority band (with 
m^/niB = 1/(1— n)) disappear at the border of magnetic- 
moment saturation rh = n. Two features are impor- 
tant here: (i) the masses are high in the almost-localized 
limit (1 — n) <^ 1, and (ii) the mass to^ in the spin- 
minority band is the heaviest, since due to the spin sub- 
band occupancy imbalance, these spin-minority quasi- 
particles scatter very strongly (due to the presence of 
the large-magnitude Hubbard term ~ U J2i 
ditional features follow from the circumstance that we 
can "switch-ofF' completely the Hubbard interaction by 
applying the magnetic field and saturating magnetically 
the system. This is possible only (and is the case) be- 
cause the field induced by the correlations hcorr enhances 
strongly the effect of applied magnetic field. In that sit- 
uation a metamagnetic transition takes also place |31]. 

The brief analysis provided above delineates the princi- 
pal message about what we mean by nonstandard quasi- 
particles in a (strongly) correlated system. First, they 
can become quite heavy, i.e. TOod/tob ^ 1. Sec- 
ond, they depend on the particle-spin direction, what 
makes them distinguishable in the quantum-mechanical 
sense, since the mass in nonrelativistic quantum mechan- 
ics is an external (input) characteristic in the problem at 
hand. Third, the effective field driven by the correla- 
tions can become very strong i.e. much stronger than 
any Weiss molecular field appearing in traditional local- 
moment magnetism. All these microscopic properties 
must be determined self-consistently. These features of 
those nonstandard quasiparticles distinguish them from 
those are defined within the original phenomenological 
Landau- Fermi-liquid theory, where the enhancement of 
the effective mass and of the magnetic susceptibility is ex- 
pressed in terms of interaction parameters. Additionally, 
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in the Landau theory of Fermi hquids the enhancement 
factors are determined by including the interaction only 
among the quasiparticles at and/or in close vicinity of the 
Fermi surface. Here, all the particles mutually influence 
each other, what is expressed via an integration over all 
occupied states when solving appropriate self-consistent 
equations. This is because in the present situation the 
interaction is strong (at least comparable) to he Fermi 
energy. 



Superconducting state: 
Fulde-Ferrell-Larkin-Ovchinnikov state 



We describe next the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state for the model of heavy- 
fermion system starting from Hamiltonian (Dl) derived 



in the Appendix D, with the hybrid pairing introduced 
in the preceding Section. 

In the standard BCS approximation with anomalous 
averages, Hamiltonian (Dl) in the narrow band limit re- 
duces to the form 



lA, 



kQ\ 



(39) 



with the single-particle energy parametrized in the tight- 
binding approximation, which in the case of square lattice 
takes the form 

Efccr = qa[—'^t{cos + COS ky) + At' cos kx cos ky],{40) 
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FIG. 9. (Color online). Phase diagram for the cases with the 
spin-dependent (a) and the spin-independent masses (b). Yel- 
low region corresponds to Q = (BCS phase), the darker one 
to Q 7^ (FFLO phase) and the white to normal state. Note 
that with increasing temperature, the transition from BCS 
to FFLO state occurs at higher fields, in qualitative agree- 
ment with experimental results. The FFLO phase is stable in 
an extended Ha-T regime only in the spin-dependent-masses 
(SDM) case i28j. 



where t and t' are the first and the second hopping inte- 
grals and is, as before, the Gutzwiller band narrowing 
factor (the inverse spin-dependent mass enhancement). 
The superconducting gap is thus determined from the 
self-consistent equation 



AfcQ = - 



N ^ 

k' 



Iklk' (/- 



k'+Qi 



fk't)- 



(41) 



Additionally, Vq is the pairing magnitude which for the 
constant hybridization has the magnitude 



Vo 



ef + U 



,1/2 



(42) 



and the factors 7^ and 7fc/ correspond to the separable 
A;-dependent factors in (D3) divided by V^. 



the matrix form 



A* 



A 



k,Q 



fkt 
k+Ql 



A2 



J2i^ki + gf^BH-^i) + N-^. (43) 
k 



The transformation to the quasiparticle representation 
has the usual form 



"fc A / Uk ""k^ ( fkt 

J {-vk ukjy 



(44) 



Next, we carry out the approximate form of the Bo- 
goliubov transformation which in the meantime acquired 
the name of the Bogoliubov-Valatin-de Gennes-Nambu 

transformation! For that purpose we represent (39) in with the Bogoliubov coherence factors given now by the 
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relations 



Vk 




'(efct + ^-k+Qi - 2/^)' + 4A2 
Cfct + g-fc+Q^ ~ 2^ 




The quasiparticle energies in the phase with Q 7^ are 



(efet + e-fc+Qi - 2/i) +4|AfcQ| 



1/2 



(47) 



where sign factor a — ± corresponds to the electron or 
hole excitations, respectively. 

Having discussed the explicit expression for the 
fermionic quasiparticle excitations, we can construct the 
free energy functional, as well as determine the system 
of self-consistent equations for AkQ,n, and m; we also 
optimize the energy with respect to the magnitude of 
the wave vector Q [33] . One should mention that for 
the electron-gas situation we have included explicitly in 
the effective Hamiltonian also the correlation field hcom 
which we optimize, whereas for the two-dimensional band 
structure ( 40 1 we have been able so far to carry out only 
the whole analysis in the Gutzwiller approximation. In 
Figs. 9 and 10 we provide the exemplary phase diagrams 
on the plane temperature T - applied magnetic field H 
for the three-dimensional gas [55] and for the square- 
lattice cases [33], respectively. Note that the phase di- 
agrams drawn in the Figures are for the simple Fulde- 
Ferrell state, i.e. for the form of the gap. Additionally 
the gap has the d-wave symmetry. 

The most important feature coming out of these Fig- 
ures are: (i) the BCS state (i.e. that with Q = 0) is quite 
robust in the lower fields and higher temperatures, (ii) 
the inclusion of the effective-mass spin dependence leads 
to a remarkable extension of regime of the FFLO stability 
in both depicted situations, (iii) in the FFLO phase the 
upper critical field of the transition to the normal phase is 
much higher than Hc2 for the BCS superconducting state 
(cf. Fig. 10, the dashed lines mark the second critical 
field Hc2 in the the Pauli limit), and (iv) the first-order 
BCS^FFLO transition can be accompanied by a weak 
metamagnetic transition [33]. One should also note that 
the detailed analysis of the FFLO state is carried out 
separately [33]. Also, the full SCA analysis of the super- 
conducting states within the full statistically-consistent 
Fukushima approach for the present model is still to be 
carried out. We do not expect though, that such analysis 
will change the picture in a decisive manner. 




0.030 



FIG. 10. (Color online). Phase boundaries for a two- 
dimensional d-wave superconductor with both the spin- 
independent masses (SIM) (a) and with the spin-dependent 
masses (SDM) (b). The FFLO-BCS transition line is of dis- 
continuous nature. The dashed line marks the stability limit 
of the BCS state as determined by the value of the second 
critical field Hc2 for the BCS state. The values of param- 
eters are n = 0.97 and Vb = 12.5 K. For these values of 
the parameters, the superconducting transition temperature 
is T = 2.5 K and the uppermost critical field for the FFLO 
phase is above 6T. Note that the FFLO state is robust in 
the situation with SDM and this result is one of the principal 
features of the present discussion [33] . 



D. Superconducting state by the Kondo-type 
(hybrid) pairing in the Anderson-Kondo model 

Above, we considered only the narrow /-band limit for 
the Anderson-Kondo model. Such model is valuable if 
the the /-level occupancy Uf can be regarded as con- 
stant in the considered regime of parameters, i.e. the 
hybridization gap is regarded as large. Here we mention. 
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FIG. 11. (Color online). Bare unrenormalized hybrid-gap 
amplitude profiles as a function of both total number of 
electrons per site or intraatomic-hybrid magnitude. The val- 
ues of parameter e/ = —1, U = 3, W — 2z\t\ = 1. The 
renormalized gap magnitude is Ac — (?(A-t-4,A4,-t-)/2. The gap 
amplitude vanishes in both Kondo-insulator {n^ — >■ 2) and in 
the localized moment {\V\ — )■ 0) limits. 



that recently we have obtained an explicit solution of 
the Anderson-Kondo Hamiltonian ( 10 1 in SCA. The de- 



tailed phase diagram will be discussed elsewhere. Here, in 
Fig. 11 we show the components A^^ = (/ifC^) and 



-A^^, both as a function of magni- 



tude V of the bare (intraatomic) hybridization and the 
total number of electrons Ue = nf + ric (per atomic site 
containing pair of orbitals / and c). The fact, that we 
treated as separate averages (f^ci) and (/^c^) and have 
obtained that they are of equal magnitude but of oppo- 
site sign means, that we have in this case indeed a pure 
spin-singlet hybrid pairing which vanishes either in the 
limit of rie = 2 (where the Kondo-insulator state be- 
comes stable) or when y — > 0, where localized-moment 
antiferromagnetic phase is stable. Also, the regime of 
AF- phase stability is separated from that, in which SC 
is stable. Work along this lines is in progress and will be 
reported separately. 



A. From real-space pairing to 
renormalized-paramagnon mediated pairing 

In this paper we overviewed the concept of real space 
pairing, induced by the kinetic exchange interaction com- 
bined with the pair correlated motion. This pairing is 
facilitated by the circumstance that the kinetic-exchange 
interaction integral J (or that for the Kondo interaction, 
J^) is comparable to the single-particle energy as ex- 
pressed by the renormalized hoping magnitude ~ \t\x 
{\V\x for hybridized systems). In that situation, the 
second exchange-coupled partner to a given electron (or 
hole) follows in their combined motion throughout the 
lattice. This situation has been animated graphically 
elsewhere [23]. To reiterate, there is no obvious interven- 
ing collective boson excitation mediating the pair bind- 
ing. A contribution of the renormalized paramagnons to 
the pairing represents an additional factor to be evalu- 
ated separately starting from RMFT. 

Namely, in the situation when the renormalized hop- 
ping magnitude can be smaller than the pairing poten- 
tial magnitude J (cf. Appendix C, when x ^ Xci), 
the Hartree-Fock-type (BCS) decoupling of the term, 
45,, « + (i3,,)4 - can be re- 

garded only as a first-order approximation. The Renor- 
malized Mean Field Theory had to be an invented as an 
approximate treatment of the effective Hamiltonians, al- 
beit not systematic in the field-theoretical sense, that we 
first perform the saddle-point approximation and include 
the quantum Gaussian fluctuations next. Those renor- 
malized Gaussian fluctuations should lead to a resid- 
ual paramagnon pairing. The last approach has been 
outlined briefly in the Appendix E. There, we reiterate 
more precisely the just mentioned difficulties connected 
with the division into the real-space and renormalized- 
paramagnon pairing parts. It remains to be seen to what 
extent the paramagnon part is important. It should be 
important in the underdoped regime x ~ Xci as the 
\t\x J. The division into the mean-field and the para- 
magnon parts may be easier to carry out in the limit of 
moderately or weak correlated systems. This is because 
then the division into the Hartree-Fock and the fluctua- 
tion parts is natural. 



B. Effective t-J and t-J-V models should include 
direct intersite (interorbital) Coulomb interactions? 



V. OUTLOOK 



Below, instead of making a summary, we pose some 
important questions concerning the real space pairing 
first and than conclude by suggesting its universal charac- 
ter applicable also to nuclear and astrophysical quantum 
matter. 



There is one additional feature of the models as rep- 
resented by the starting effective Hamiltonians ^ and 
([5]) not discussed in detail. Namely, the respective in- 
tersite or interorbital direct Coulomb interaction should 
be included in the corresponding Hamiltonians. This is 
because, from one side J < i^(ij) and < Ufc, and 
from the other, such terms are already included in the 
corresponding full expression of the Dirac spin-exchange 
operators. Their role many not be so obvious if we de- 
couple them and include anomalous contraction of the 
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type (aicraia)- One may think, that when decouphng in 
the Hartree-Fock manner the term ~ PiVj (or vinm in 
the hybrid pairing case), the contribution to the singlet- 
and the triplet-pairing channels is the same for this spin- 
independent term and therefore, the whole term can be 
disregarded when considering the paired state. Addi- 
tional reason f or disregarding the term ~ K^j is that the 
paring is due to the combined effect of exchange and cor- 
related hopping. Nonetheless, the inclusion of such terms 
would present itself as an additional test of the whole pic- 
ture, though most probably the expected results would 
be of secondary importance in the most important sit- 
uations (c.f. also the related discussion of the situation 
with Kij < in the main text). 
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C. Prom real-space pairing to the Hund's rule 
exchange pairing 

Once we have suggested, that the exchange interac- 
tion is not only a source of spin magnetism but also a 
fundamental mechanism of the real-space superconduct- 
ing pairing, we may ask if other exchange interactions 
may become a source of pairing or superfluidity. Here, we 
have in view our ideas [31] about the role of the ferromag- 
netic, intraatomic, interorbital (Hund's rule) exchange in 
orbitally degenerate systems. Customarily, the Hund's 
rule exchange is associated with itinerant ferromagnetism 
in moderately correlated systems [35, or with an orbital 
ordering mixed with magnetic ordering in strongly cor- 
related systems [2]. We have suggested that for the itin- 
erant correlated and orbitally degenerate narrow band 
systems the Hund's rule exchange provides a stable spin- 
triplet superconducting state both in the Hartree-Fock 
[55] and in the strong-correlation limits [37] ■ Fig. 
12 we illustrate this statement in the former limit, re- 
spectively by plotting the corresponding phase diagram. 
The phase diagram displayed here supplements the well 
known phase diagrams involving only magnetic phases 
[36] with the superconducting paired states. This topics 
require a detailed separate analysis, which is in progress. 



FIG. 12. (Color online) Exemplary phase diagram [35] on the 
plane band filling n - Hund's rule exchange J = . The 
phases are as follows: NS - normal state; A -superconducting 
spin-triplet state with the gaps A^^ = An; Ai - supercon- 
ducting spin-triplet state with the gap A^^ 7^ only; AF 
- antiferromagnetic state; FM - ferromagnetic state; SFM - 
saturated ferromagnetic state; AF+SC - mixed antiferromag- 
netic -I- superconducting state with the superconducting gaps 
different for majority- and minority-spin electrons on given 
site. The square lattice was assumed with W = 1, U = 7J, 
and no hybridization in the doubly degenerate band was in- 
cluded. 

dressed property within RMFT, i.e. without the phase 
fluctuations of the superconducting gap. 

One additional general remark is in place at the end. 
Namely, the exchange interaction is a universal interac- 
tion in the sense, that it takes place between any two 
interacting fermions. Therefore, it would be interesting 
to address the question of real space pairing in a corre- 
lated quantum nuclear or astrophysical matter (nucleons 
in nuclei, neutron stars, quark-gluon plasma, etc.). This 
topic could become a further example of incorporating 
the laboratory condensed matter physics into fundamen- 
tal quantum-matter physics. 



D. Concluding remarks 

The real space pairing concept is different from other 
concepts, which are based on the idea of virtual boson 
excitation (phonon, paramagnon etc.) mediating the at- 
traction between fermion. Here the exchange interaction 
combined with a correlated motion of the pair (pair hop- 
ping) is responsible for the pairing. The latter idea has 
an intuitive interpretation (23| . However, the real space 
pairing concepts should be tested further. A selected 
comparison of RMFT results with experiment have been 
also discussed in this text. Those results are insufficient, 
as for instance the Fermi surface evolution and the pseu- 
dogap appearance have not been (and cannot be) ad- 
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Appendix A: The concept of t-J model and 
real-space pairing - in perspective 

In 2011 we celebrated 25-th anniversary of the discov- 
ery of high-T(7 superconductivity. In 2012 there is the 
35-th anniversary of our publication of t-J model [51] , as 
well as 25-th anniversary of discovery its meaning in the 
context of high-temperature superconductivity. Below I 
provide a somewhat biased personal account of the last 
two topics. 



1976-1986 

The generalization of the ideas of Anderson [31 [35] con- 
cerning the origin of antiferromagnetic (kinetic) exchange 
interaction for the Mott insulators to the strongly corre- 
lated metals was was done by the present author in 1976 
and subsequently published with his colleagues [39] . The 
pioneering period was reviewed in my habilitation-schrift 
[40] . The idea of defining the division of the hoppings 
into the four terms corresponding to the intra- and inter- 
Hubbard subband hopping was particulary difficult to 
envisage. So the resulting canonical transformation, as 
the "unperturbed" part of the Hamiltonian, was of non- 
diagonal form in the Fock space. The last feature meant 
that we had to perform the canonical perturbation ex- 
pansion in the operator form. A more formal overview 
of the pioneering era is given elsewhere j41|. Many au- 
thors have rederived subsequently of the t-J model, but I 
claim that our derivation was the first one, albeit limited 
to the discussion of magnetic phases and the Mott tran- 
sition as the band filling approaches unity. Some authors 
claim that Harris and Lange [12] provided the basis of 
projected fermion operators {ajg.(l — nia),aia{^ — nio-)} 
and the kinetic exchange for the n ^ 1 case, but a cur- 
sory look to the paper shows that it does not contain 
any explicit mention of the kinetic exchange and does 
not introduce t^ /U effects explicitly. 

One should mention that practically at the same 
time we have also introduced [13] an analogical effective 
Hamiltonian for the Wolff model of the magnetic impu- 
rity, which is to a certain degree analogous the Anderson 
model of magnetic impurity. There, we introduced the 
cases of " shallow" and " deep" impurity cases, which lead 
me to the concept of the modified Schrieffer- Wolff trans- 
formation in the "shallow-impurity" case. This distinc- 
tion is elaborated in the main text, when we talk about 
the hybrid (Kondo) pairing for the Anderson lattice, ap- 
pearing concomitantly with the itineracy of / electrons. 

Practically, nobody was interested in those papers 
then, as it seemed that the fashionable mixed-valence 



and heavy-fermion physics had not much to do with the 
Hubbard model or the Wolff model. There, the peri- 
odic Anderson model or the Kondo-lattice models were 
regarded as the distinct and relevant models. 

1987-now 

The revolution in theory of correlated electron sys- 
tems came with the introduction of real-space pairing 
amplitudes {a\^a\g) within the t-J model. The author 
learnt about the idea from the preprint of Ruckenstein, 
Hirschfeld, and Appel [H^- The often quoted paper [IS] 
was completely illegible at the time (not only to the au- 
thor) and it is extremely difficult now to argue, at least 
for me, who introduced the real space pairing first. We 
can also say, that we have invented [B] a correct form of 
the t-J Hamiltonian with the precise real-space projected 
operators {Bjj, Bij}, as well as have extended the repre- 
sentation to the case with projected hybrid pairing op- 
erators {bljj^,bim}, which appear in the Anderson-lattice 
model in " large-but-finite U limit" [7] . The last model we 
have termed as the Anderson-Kondo model. But, in anal- 
ogy to t-J model, this model should be rather termed the 
t-J- V model, where V stands for hybridization which ap- 
pears concomitantly with the exchange couplings (Kondo 
or superexchange) and is instrumental in driving the itin- 
eracy of / electrons in heavy-fermion systems. 

The second aspect of the real space pairing driven 
by the kinetic exchange interaction is the absence of 
the virtual boson driving the real space at least in the 
mean field approximation. This question leads us to a 
highly nontrivial problem of formulating renormalized 
mean field theory (RMFT). This nontriviality of the 
mean-ficld-approximation formulation in this situation 
stems from the fact that the projected hopping part is, 
strictly speaking, of many-body nature, since it contains 
al^{l — nig)ajcr{l — nja) factor. However, absolutely 
crucial in the analysis are also the statistically consis- 
tency conditions phrased explicitly in our recent works 
[5] [T31 [TJ, as well as in some other papers [15]. The 
rationale behind these consistency conditions, expressed 
through Lagrange multipliers and added to the effective 
Hamiltonian, is to ensure that the self-consistent equa- 
tions for the averages appearing in the RMFT p rovide 
the same results as those obtained from an appropriate 
variational procedure for the free-energy functional rep- 
resentation, the Landau functional in this situation. 

The third aspect of the current research is to go beyond 
the mean-field approximation, i.e. include the Gaussian 
fiuctuations. So far, this approach has been formulated 
within the slave-boson approach (SBA) [3], which how- 
ever contains spurious phase transitions corresponding to 
the condensation of auxiliary boson fields. RMFT with 
the constraints is in some aspects equivalent to SBA in 
the saddle-point approximation. A systematic approach 
to incorporate the quantum fluctuations starting from 
RMFT, is still missing [cf. Appendix E]. 
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An approach based on the Quantum Monte-Carlo 
method provides a very important insights into the re- 
sults, albeit limited to very small systems. A combina- 
tion of RMFT with spatially inhomogenous order param- 
eters and quantum Monte-Carlo methods seem to be also 
very promising. 



Appendix B: Analytic estimate of the BCS type gap 
magnitude, universal 2A(0)/fcsTc ratio, and 
reduction to classic mean field case 



Here we sketch the derivation of Eq. (|32j). Its im- 
portance derives from the surprisingly good approach in 
Fig. 6. of the numerically obtained gap magnitude Ac 
with that from (32 1. Our results put on a solid ground 



the estimates proposed originally (but not derived) by 
Rickayzen PSI . 

We start from the self-consistent equation for the gap 
magnitude within BCS theory, namely 



1 - 



— hcor 



dep{e)V{e) 



1 



Ve2 + A2 



: tanh 



Ve^TA2 

2knT 



(Bl) 

where ^(e) is the absolute value of the pairing potential, 
p{e) is the density of states in the band, e is the parti- 
cle energy counted from the Fermi energy e^? = 0, and 
is the energy cut-off for the pairing potential (in the 
real-space pairing case, we take hwr, ~ J, the magnitude 
of kinetic-exchange interaction). We believe, it is rea- 
sonable to make the so-called BCS approximation at this 
stage, justified in the present situation in the following 
manner. Namely, we utilize RoUes's theorem about the 
average value of the defined integral: 



dx f{x) = f{x){b-a), 



(B2) 



where x € [a, 6]. In other words, the area under the 
curve f{x) can be represented by that of a rectangular. 
In applying this theorem to (Bl) we find that it can be 
rewritten in the form: 



1 = hLUDPi^) 



1 



tanh 



x/e2 + A2 
2k bT 



(B3) 



Now, it is reasonable to assume that A^/ia;£)<CeF,so 
we can represent — /Lt in the same manner, as in the 
Fermi liquid theory, i.e., assume that near the Fermi level 



Ek - p = hvpk, 
so that in volume ft the density of states is 

n 



with e — hvpk, 



(B4) 



(B5) 



(hvp) 7r2 

which is an antisymmetric function of e. Hence, the in- 



tegral of ( B2 1 is an increasing function as we cross the 



Fermi level energy p. This minor argument is used to 
make a bold statement, that in such a situation an al- 
most exact approximation can be made , i.e. one can 
put e = in ( B3 ) . This means that the integral equation 



can be replaced by the following algebraic equation 



1 = hujDViO)p{0) 



1 



A(r) 



tanh 



2kBT 



(B6) 



which represents a classical mean-field self-consistent 
equation for the order parameter A(r) (e.g. for the 
spin (1/2) Ising model). Parenthetically, this means that 
the BCS-superconductor phase transition is of classical 
nature, with a mean-field (van-der-Waals type) classical 
critical point Tq- 

One can draw additional conclusions from the result 
(B6). First, assuming that the order parameter A — 
A(r) is small when T -> Tc — 0, one can rewrite ( B6 1 in 



the form of Landau equation for the order parameter: 



A ~ hwDV{0)p{0) 



or equivalent (for T ~ Tc): 



A 



1 



3 \2kBT 



2kBT 



A 



hujBVp 

2kBT 



1 fvjj£,V p 
24 {kBTcf 



A2 =0. 



Thus, to a good approximation 

kBTc = ]^iODV{Q)p{Q) 
and then it can be rewritten in the form 



A(r) 

2kBTc 



= tanh 



A(T) ' 
2kBT 



(B7) 

(B8) 

(B9) 
(BIO) 



This relation, in turn, leads to the following universal 
relation for the gap-to-Tc ratio: 



2A(0) 



4. 



(Bll) 



kBTc 

This value differs from the usual estimate of this ratio 
- the value 3.53 (or 3.311) - obtained directly from Eq. 
(Bl ) estimate. We think the discrepancy is not large , in 



view of even larger difference in experimental data for the 
BCS superconductors. In the context of our present dis- 
cussion of correlated systems, this ratio of maximal am- 
plitude At,(0) to the real experimental value of Tc should 
be either equal to the (Bll) or even larger. The latter 



situation takes place if the pairing-potential renormal- 
ization is augmented by the additional renormalization 
of the density of states, so that 2A(0)/Tc = 4(g*)~^ 
On the basic of the above treatment one can also 



rewrite equation (B6| in the form 



A(T) 
A(0) 



tanh 



A(T) Tc 
A(0) T 



(B12) 



which coincides with that proposed originally by Rick- 
ayzen |26j . In brief, this discussion, ignited by the coin- 
cidence of the results shown in Fig. 6, illuminates a "BCS 
aspect" of the present renormalized mean field theory. 
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Appendix C: Elementary estimate of the lower 
critical concentration for onset of 
antiferromagnetism 

Here we estimate the critical concentration for the an- 
tiferromagnetism disappearance as a function of doping. 
For that purpose we start from the canonical version of 
the t-J model ^ without the last two terms included. 
One may say that in the Mott insulating state antifer- 
romagnetic interaction dominates over the hole hopping 
and this means that the (lower) critical concentration 
is achieved when the two contributions are of the same 
magnitude, i.e., 



z\t\x{l-x) + J{l-xf 



0. 



This happens roughly at the hole concentration 



J 



1 J 



z\t\ +J z \t\ 



(CI) 



(C2) 



which yields ~ 0.07, a quite reasonable value in view 
of the simplicity of the estimation. 

In order to try to improve the estimate we formulate 
an elementary Gutzwiller type of approach in which the 
renormalizcd hopping and the exchange parts are of the 
same amplitude. This means that we can postulate the 
ground energy per site of the form 



Eg 
N 



= ~z \t\ $(A) V(aLaj<x)o - JzX 



(C3) 



where the variational parameter A = {{l/4:)ninj—Si-Sj), 
is equal to unity in the Neel state of the Mott insulator 
and reduces to the value A = {l/A){ni){nj) = 1/4 in the 
uncorrelated state. The renormalization factor can be 
expanded in powers of A 



= .90 +5i'^ + ff2A^, 



(C4) 



Hence, 



1 



8|t| 
3J 



0.1, 



(C8) 



for \t\/J — 3. Finally, the optional ground state energy 
per pair of states and for n ^ rid i*^ 



— ^-\t\n{l-n)~-n 



64 \t\ (1 - n) 



(C9) 



Note that if the ratio \t\/J = 6, then Xd ^ 0.07. This 
is still too high value by a factor of 2 when compared to 
the experimental value: Xc-^ ~ 0.03 for La2-xSr^Cu04. 



Nevertheless, the expression (C8) has a straightforward 
interpretation. Namely, the first term is the renormal- 
ized band energy in the Gutzwiller approximation in the 
normal state and in the U = oo (J = 0) limit. The sec- 
ond term is the energy of the Neel state in the mean field 
approximation for the Heisenberg part. The last contri- 
bution is the higher-order contribution coming from the 
correlations (note that our order parameter A in the Lan- 
dau expansion ( C3 ) is of two-particle nature, but still we 



recover a single-particle (mean- field) description). 

This lower concentration (or Xc-i ) is missing in the 
present approach (cf. Sec. IV C), since no localiza- 
tion effects within RMFT-SCA were taken into consid- 
eration. This means that we may need to supplement 
our whole approach with some sort of cluster expansion, 
i.e. improve on the Gutzwiller (or Fukushima) ansatz. 
We should see a progress along this line in near feature 
e.g. by evaluating the averages with the help of a full 
Gutzwiller wave function 1221. 



Appendix D: From Kondo coupling to real space 
pairing in an almost localized Fermi liquid of heavy 
quasiparticles 



in which higher order terms have been disregarded as 
we are close to the magnetism disappearance (or onset). 
Also, the hopping probability in the uncorrelated state 
is {al^aja)o = ncr(l - n^) = ("■/2)(1 - n/2). Next, the 
constant go, gi, and 52 can be determined from the par- 
ticular solvable limits, namely: 

1° $(A = 1) = (C5) 
we do not encounter any hopping, i.e., (?o + 5i + 32 = 
2° $(A =l) = l (C6) 



i.e., 50 + hi + T692 = 1- 



Hence, the optimal value of A = 1/4 (from dEo/dX = 0) 
becomes 



A 



1 



Jn 



4ff2 \t\ (1 - n) 



91- 



(C7) 



Now, the remaining condition is that A = 1/4 when J = 
0. In result, we obtain go = —-^t gi — ^, and 52 = ~if • 



Hamiltonian ( 10 ) can be transformed to the mean-field 



form and subsequently to the momentum representation, 
in which it takes the form 1331 



ka 

f kk'Q 



(Dl) 



with 



^kQ — ^2 \-'k+Q/2t-k+Q/2l 



Jk+Q/2l^~k+Q/2t )■ 



(D2) 

The wave vector Q is nonzero in FFLO phase. The single- 
particle part can be easily diagonalized by moving to the 
hybridized basis, whereas the pairing part is represented 
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by a separable pairing potential. We shall proceed with 
the transformation to the hybridized basis first, which 
yields the following transformed pairing part for the lower 
hybridized band 



ef + U 



E 

fcfc'Q 



\VkVk'\H<iaq-, 

^ "fc+Q/2t"-fe+Q/2i"-fe'+Q/2i"fc' + Q/2T' (D3) 

with the dispersion relation in the lower hybridized band 
in the form 



Eka 



(D4) 



Note that in (D3| and (D4) we have written the formu- 



las for the spin-polarized situation. For real Vfe, the hy- 
bridized quasiparticle operator reads aka — cos Oka fka + 
sin6'feo-Cfccr) with the condition for the mixing angle 9k 



tuation parts according to 

?i«HGA-A(")-^(5,-(5,)) 

i 

i(T 

-^(a.-A(™)).(5.-(5.)) 

i 
i 

where we have assumed that {Si) is oriented along z- 
axis, if A(") = (0, 0, A^™)) is taken. The first three terms 
compose the SGA Hamiltonian ( 18 ) (the factor 1/2 in the 
second term is irrelevant, since we can redefine A'™-* = 
I a'™-*), whereas the remaining two terms represent the 
fluctuation around the mean-field equilibrium. We can 
write down the fluctuation part in the form 



(E4) 



tan 20ka 



Wka 



(D5) 



Taking the states on the Fermi surface we have 
ta,n26ka ~ 29kcr, and hence 9kcr = Vka/^fa- In effect, 
aka — fka- The last estimate provides us with the 
starting Hamiltonian (39 1 in the single-band limit. The 



complicated fe-dependence of the pairing potential means 
that in general, the nature of the superconducting gap 
may take a form more complicated than pure extended 
s-wave or d-wave forms. Note also that even though the 
pairing potential contains explicitly the spin quantum 
numbers, it is in fact spin independent, as assumed in 
the main text. 



Appendix E: Incorporation of quantum fluctuations 
within the SGA - RMFT approach and the eff'ective 
fermion-boson model 



The effective Hamiltonian ( 18 1 in SGA approximation 



can be generalized to include the quantum fluctuations. 
Namely, we can generalize it the spin rotationally invari- 
ant form, which can be now written in the following real- 
space form: 

n « UGA-Y^Xt"^ ■ (5, - (5,)) - ^ Ai") (n,, - (n,,)) , 

(El) 



where 



= \ Y'y'^),yo' alalia' 



(E2) 



and r = (r^,, r, r^) are the Pauli matrices. The con- 
straints can be divided into the mean-field and the flue- 



It is composed of the local spin and the charge fluctu- 
ations coupled to the fluctuating Bose fields 5A|™' and 

(5Af \ Note that we have neglected the fluctuations in 
double-occupancy df = {ni^rii^), which can also be in- 
cluded, at least in principle. The form (E4) represents 



the starting point of analysis, which will not be detailed 
here. 



From t-J to the fermion-boson model with 
correlations: A formal transformation 

The above intuitive picture can be rephrased in a 
more formal language. We outline this method for the 
t-J model containing explicitly the real space pairing. 
Namely, we start from the expression for the partition 
bmction of t-J model (for the simplest situation, i.e., with 
no three-site terms and with Kij = 0) in the coherent- 
state representation, which reads 



Z = D 



exp \- f drY, [4(r) [{d/dr ~ m) % + U,] 

6,,(r) + ff,(T)]|, (E5) 

where we integrate over the Grassmann fields {b\^,hjc), 
P = {kBT)~^ is the inverse temperature, and 



(n) 



(E6) 
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is the interaction part. Now, the novel basic idea is that 
the summation over bonds (ij) can be regarded formally 
as a simple summation; B\-Bij then has the form ap- 
propriate for the Hubbard-Stratonovich transformation, 
i.e., we can linearize the quartic term in the following 
manner: 



(ij) ^ 



A- -A*- ' 

(E7) 

9 



with g = J(^ij^ — J. In effect, the partition function takes 
the form 



IA-12 



E 

(iJ) 



A,,(r)4(r) + A*,(r)i?,,(r)-^ 



(E8) 



In effect, we have the following Hamiltonian when we 
include both fields: 

ija 

-5][A,,(r)4(r) + A*,i3,,(7 



^JriA^(r)A,,(r). 



(E9) 



We see that Ay (r) in the RMFT is played by A,j = 
Jij{Bij). Here, the representation is exact and the "time 
r" dependent pair-fermion field Blj{T) is coupled directly 
to the bosonic field Aij(r) [17], the latter is a Gaussian 
fluctuating complex field. Alternatively, one may add 



^^j^^ A*jAij to the effective Hamil- 



the Bose term — J2{ 

tonian as in ( |E9[ ). In the saddle-point approximation the 
T dependence of the fields is disregarded and the pairing 
amplitudes are minimized for the free-energy functional. 
Note, that A^ represent the spin-singlet-pairing ampli- 
tudes, as bIj and Bij have the form mi. However, even 
in the mean-field approximation, defined in the above 
manner, the problem is still highly nontrivial, since the 
first term in (|E9|) contains the projected fermionic opera- 



tors bff^'^^ and b^rr with non-fermionic anticommutation 
relations ([3]) . This was actually the reason for the discus- 
sion on nontriviality of the renormalized mean-field the- 
ory (RMFT) . This circumstance leads also to the conclu- 
sion, that if we start our analysis from RMFT for ( E9 1 , 



then to be able to include the fluctuations in Aij in a 
systematic manner, one has to include the fluctuations 
in projected fermionic fields around RMFT solution as 
well, which may represent a formidable task. Most prob- 
ably, the viable method to diagonalize exactly (E9| in 



the saddle-point approximation is to resort to the Quan- 
tum Monte-Carlo or related cluster-expansion methods. 



Nonetheless, the expression (E8) provides a hint concern 



ing the direction, in which a systematic approach towards 
rigorous solution of the t-J model should be developed. 
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